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Abs  t  rae  t 


(.fleets  of  Laf:ii'tic  shear  on  lower  hybrid  inodes  are  investigated.  It 
shown  that  due  to  r.on-lt  cal  effects,  even  a  snail  slit  nr  can  si.vii  f  i  cant  ly 
affect  tlie  instability,  leading  to  stabilization  for  sou  par. ureter  r <ing.es 
These  results  are  of  importance  in  the  context  of  the  recent !y  proposed 
mechanism  of  lower  hybrid  acceleration  and  ion  evolution  '■  the  swrraurora 
r  e  g  i  o  n . 


I.  Introduction. 


It  has  recently  been  shown^  that  ions  can  be  accelerated  perpendicularly 
to  the  magnetic  field  by  resonant  interactions  with  current  driven  lower 
hybrid  modes.  This  accelerated  portion  of  the  ions  can  evolve  into  conic 
distributions  and  propagate  upwards  along  magnetic  field  lines.  When  they 
reach  the  region  where  they  can  be  strongly  energized  by  the  electrostatic 
shucks,  it  is  argued^  that  the  resulting  ion  distribution  can  lead  to  the 
excitation  of  electrostatic  ion  cyclotron  modes.  This  provides  a  plausible 
explanation  of  the  simultaneous  observation  of  the  electrostatic  ion  cyclotron 
modes  with  the  kev  ion  distributions  in  the  suprauroral  region. 

Since  the  driving  current  also  produces  a  magnetic  shear,  which  generally 
exerts  a  stabilizing  influence,  it  is  quite  important  to  investigate  whether 
the  lower  hybrid  mode  remains  unstable  inspite  of  the  influence  of  shear.  The 
magnitude  of  the  shear  is  generally  quite  small,  with  shear  length  Ls  of  the 
order  of  5 i .1  ( >  km,  and  it  may  seem  reasonable  to  ignore  it  altogether.  However, 
our  recent  studies-*^  in  the  context  of  the  current  driven  ion  cyclotron 
instability  have  shown  that  even  a  small  shear  can,  due  to  non-local  effects, 
sometimes  produce  a  very  significant  reduction  in  the  growth  rate  of  an 
instability;  under  some  conditions  it  even  leads  to  stability. 

Such  an  investigation  was  begun  during  the  summer  program  and  led  to 
interesting  preliminary  results.^  This  effort  was  continued  under  the 
milligram  program  and  the  results  todate  indicate  that  the  effect  of  magnetic 
shear  can  be  very  significant,  and  can  lead  to  a  stabilization  of  the  lower 
hybrid  mode  fur  a  wide  range  of  physical  parameters  ol  interest. 


Si  i  e  ri  L  i  f  ic  Backprou  n d  . 


11  . 


Tht1  Importance  of  lower  hybrid  wave*  excited  by  energetic  electron  beans 
in  tb.e  context  of  ion  acceleration  processes  has  recently  been  pointed  out,^ 
and  it  has  been  argued  that  a  transfer  of  energy  from  the  electrons  to  the 
ions  is  effected  due  to  the  simultaneous  resonance  of  the  lower  hybrid  (LH) 
mode  with  both  the  electron  and  ion  populations.  Tins  ion  acceleration 
process  is  particularly  efficient  at  lower  altitudes  where  the  lower  hybrid 
waves  have  a  high  intensity  over  a  broadband  of  wavelengths.  Typically,  1  eV 
ions  can  bo  raised  in  energy  to  hundreds  of  eV  or  beyond  in  the  suprauroral 
region  by  the  current-driven  LH  modes.  Since  the  ions  are  energized  primarily 
transverse  to  the  field  lines,  they  acquire  "conic”  distributions,  with  pitch 
angles  clustering  around  90°  to  140°.  Because  of  the  mirror  geometry  of  the 
earth's  magnetic  field,  the  transverse  energy  gained  by  the  ions  will  be 
converted  to  longitudinal  energy  as  they  move  upwards.  When  the  ions 
propagate  across  a  kilovolt  electrostatic  shock,  they  become  field  aligned. 

The  combination  of  such  keV  ion  distributions  and  of  that  of  the  background 
ions  can  lead  to  the  excitation  of  electrostatic  ion  cyclotron  modes.  Thus, 
t he  LI:  acceleration  mechanism  may  provide  a  possible  explanation  of  the 
simultaneous  observation  of  the  electrostatic  ion  cyclotron  modes  with  the  keV 
ion  distributions  in  the  suprauroral  region,  particularv  at  altitudes  above 
i , '  on  km. 

Central  to  this  scenario  is  the  idea  that  high  intensity  broad  band  LH 
waves  are  excited  by  the  energetic  electron  hems.  However,  it  should  be  noted 
that  this  conclusion  is  based  on  the  simplified  electrostatic  dispersion 
re  la  t  i"ri ,  1  ^  and  the  effects  of  the'  mngne'tic  shear  generated  by  the  field 

ali  gr.e-d  currents  have  not  been  considered  in  equation  (2)  of  Reference  1  or 


equation  (1)  of  keierence  5.  Magnetic  shear  general ly  exerts  a  stabilizing, 
influence,  and  if  it  could  effectively  reduce  the  growth  rates  of  the  LU 
tr  odes ,  the  above  mentioned  explanations  of  the  ion  acceleration  process  and 
the  generation  of  electrostatic  ion  cyclotron  modes  would  become  untenable. 

Thus  recognizing  the  importance  of  investigating  the  influence  of  shear  on  the 
current-driven  Lit  modes ,  we  have  initiated  a  systematic  study  of  tnis  problem. 

We  were  motivated  to  examine  this  question  inspite  of  the  small  strength 
of  shear  because  of  past  experience  and  background  in  the  context  of  our 

0  Q 

extensive  study->J  of  the  current  driven  ion  cyclotron  instability,  where  we 
have  shown  that  even  a  small  shear  can,  due  to  non-local  effects,  produce  a 
significant  reduction  in  the  grow’th  rate  and  also  reduce  the  band  width  of  the 
wave  numbers  for  which  the  mode  is  unstable. 

The  formulation  of  the  problem  of  the  lower  hybrid  mode  in  the  presence  of 
magnetic  shear  is  given  in  the  next  section,  followed  by  the  method  of  solution 
and  <;  discussion  of  results  in  subsequent  sections. 


111.  Korn:  1 .1 1  i  wn  of  tin.  Shear  I'rohler. 


For  i  li-ct  ros  tat  i  r  waves,  the  dispersion  relation  for  tin  current-driven 
lower  hybrid  mode  is  given  by  1 ^ 
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electron  density,  Vtb  indicates  the  thermal  spread  of  the  beam  around  the  bear. 

2  2 

velocity  Lt,b,  and  we  have  assumed  Wp0  <<  wce. 

The  solution  oi  the  1 o c a  1  dispersion  relation  (1)  is  given  implicitly  by 
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Fo i  (f/u)  <<  1,  the  local  dispersion  relation  sinpli  f  it-s  t< 
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The  non-local  formulation  Is  achieved-*^  by  introducting  (i)  k  =  kv(x/Ls), 
where  x  is  tilt  distance  measured  perpendicular  to  the  current  sheet  and  ]  is 


the  shear  length  and  (ii)  kx  = 


,  which  converts  F.q .  (1)  into  a 


differential  equation.  The  magnetic  field  is  along  the  z-direction  and  k  = 

1 

?  2  2  dz 

kv  +  k^  -*■  ky  -  ^—7 — .  The  resulting  differential  equation  alone  with 
appropriate  boundary  conditions  constitutes  an  eigenvalue  problem  for  the 
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With  these  substitutions,  Eq .  (1)  can  be  rewritten  as 


,  ri  .2  d‘  ,  ,  .  , 

<TT)  — 


where,  invoking  !d  <<  1,  we  have 
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Tiie  local  dispersion  relation  is  recovered  from  tlq.  (7)  by  taking  Ls 


i.e. 
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which  is  easily  verified  to  lead  to  hqs.  (2)  and  (5). 


i  V .  Solution  o*  tlit 


Shear  Prntla. 


'it.'  (i<  y,,rturt-  f  roc!  the  local  dlspt-rti  ion  relation  0  =  o  cat;  be  invest  igateri 

Kv  m  •  v  1 1. ••  tie.  differential  equation,  (7),  subject  to  the  boundary  conditions 
t  hut  the  electrostatic  potential  function  regains  bounded,  or  has  outgo! ng 
etieryv  b>;::;.larv  condition,.  In  gtntral  the  solution  to  the  differential  equation 
cannot  K  obtain'd  in  dosed  fort.,  and  numerical  solution  or  a  computer  nay  bi 
required.  However,  one  cast  take  advantage  of  the  physical  characteristics 
embodied  in  the  i>  function,  and  in  particular,  use  a  local  Taylor  expansion 
a  round  sore  ancle  ,  to  be  specified  later,  to  obtain  an  approximate  form  for  n 
which  enables  cue  to  obtain  an  analytical  solution.  This  is  the  approach 
adopted  litre. 

out  can  rewrite  Kq .  (7)  in  the  form 
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for  small  (a. la),  ')  is  essentially  a  parabola  and  this  suggests  expanding  0  upto 
second  order  ar  >u:id  some  convenient,  is  yet  arbitrary  "angle"  v0, 
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i ; .  e  i..; .  u;  j  «  a  We  her  equal  ion,  with:  a  set  of  solutions  given  by  the  product  of 
a  C'oissi.ri  and,  n’erri  te  polynomials  and  an  eigenvalue  equation  tor  w 
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i  re.qui  rener.t  s  of  self-consistency  of  this  approach  is  that  the  cubic 
’  urdi r  terms  not  included  in  Eq .  (11)  he  smaller  than  the  terms  that 
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he  weak  shear  liir.it,  Lg  »>  r  j ,  one  can  use  the  sinpli  r  Eqs.  (13)  and 
exp  licit  evaluation  of  o'  and  based  on  Eq .  ( l(')  leads  to 
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t>  r  :•  =  U‘,(.  k‘  Vlb)  =  (i  j /T(  y  k'  r\  ,  \  =  '*ei  ,  nay  rar.Kt- 

.  .j-'-eri  V  the  results  f  -r  i  =  In.  for  the  central  r  arty  os  of 
,  i.  -  1  t  •  •  t  or  :  ~  : 1 . 1 1  i  a:id  !-.l  r>  ‘••pent  i  ve  !y . 

1  :  :..•  characteristics: 

■  ’■  ;•  ,  t  here  is  < .  I  v  me  ramie,  which  is  da~-pi  d.  T*  < 

s  :  •  ■  <■  :  r  •  ■ :  r  ;  1  (  i  •  i' .  w  i  t 1 1  s  r;a  1  1.  e  r  a  )  . 

.  .  i  •  .  :  ■ ,  r  -  :  :  i:  .;:  I  st  rone  hear  st  ron  ft h  (i:  1 1  ■  < 

•  *.  •  It  :t  :  1  i-'.rm  :  >:  t * •  <  Table).  7:.<  ’”  o’nitnd>  of  • 

.'.I  a  .  t  •  r-  rive  a  yrowinc  r  -c!<  (e.y.  u  = 

'  .  i  ' . : ;  i  t  ' .  ■  -.  i  .  ;  alt  is  t  1 ;  t  t!  t  S',  i  1 1. 1 : .  t  O:  n  W!.i  :  t  it 


t  ■  r- 


r*-i  at:!  lead-  to  a  d  nan  n  w!i»-ri  no  nor"::, 
ar  sr .allot).  Cor  respond i ;  '  rami",  uncut 
1 1 1  r.ot  covered  tv  the  r.m.'e  y  i  von  i ::  tin 
=  •  * ,  i  =  k.  3  i ;  'a  =  3 ,  '.••  ~  s ;  l  =  3 ,  i  - 
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Table  1 

(5  =  0.01 

Y  (In  lu“4) 

5 

4 

3 

-0.  57 

-0.  72 

-0.72 

-0.6  5 

-0.83 

-0.  94 

-0.76 

-0.99 

-1.25 

-1.02 

-1.28 

-1.70 

-2.05 

-2.23 

-2.73 

+0. 01b 

-2.58 

-2.69 

-3. 15 

-0.  15 

+0.072 

-3.75 

-3.73 

-4.06 

-1.18 

+U.022 

XXX 

-4.80 

-5.05 

-0.58 

+0.063 

(A )  For  very  large  u,  (e.g.  u  =  10) ,  a  fourth  dona  in  appears  for  strong 
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bear.,  strengths  (a  <  10).  This  is  a  growing  mode  with  a  significant  >  ,  somewhat 
below-  the  local  growth  rate.  However  this  domain  is  too  far  away  from  the 
parameter  region  of  practical  interest.  By  extrapolation  such  a  region  might 
occur  for  moderate  u  values  if  we  make  a  sufficiently  small,  below  a  =  1.  But 
a  =  1  corresponds  to  the  case  where  the  beam  density  is  the  same  as  the  ambient 
electron  density, 

(5)  For  large  a,  (i.e.  weak  beam  strengths)  one  can  show  analytically  that 
:  scales  as  a-*-.  This  has  been  observed  in  the  results  for  large  a  ~  1(1^.  The 
a.-y;:jtotic  scaling  sets  in  for  smaller  a  values  for  smaller  u. 

to)  i  is  proportional  to  fc  in  this  range,  as  can  be  seen  by  comparing,  Tables 
1  and  2. 

From  these  observations  we  conclude  that  in  the  region  of  practical 
interest,  there  is  only  the  damped  mode  which  corresponds  to  the  damped  branch  of 
t  :.e  linear  theory.  The  second  mode  (which  is  marginal ly  positive  or  negative) 
occurs  only  for  strong  currents.  Thus  thi  current  driven  lower  hybrid  mode  is 
essentially  stabilized. 
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Ue  haw  formulated  the  problem  01  the  current  driven  1  owe r  hybrid  n:«  Jt  in 
tin  presence  of  shear  and  developed  the  ntce-ssary  m.  a-local  t  rt  atm,  nt  for  its 
solution.  Within  the  assumptions  of  the  problem,  wt  find  that  the  mode  is 
t  sseat i ally  stabilized.  This  is  a  si tnif icant  result,  especially  in  the  context 
or  the  proposed  ion  acceleration  mechanisms  based  on  the  lower  hybrid 
instability.  In  view  of  the  importance  of  the  result,  we  comment  briefly  on  the 
assumptions  and  the  mode  of  solution,  ami  indicate  further  avenues  of  extending 
t!iis  work. 

We  have  started  with  a  simplified,  electrostatic  local  dispersion  relation. 
The  full  electrostatic  dispersion  relation  expressed  as  an  infinite  sum  over  the 
ion  cyclotron  harmonic  terms  should  be  employed  as  the  starting  point. 

The  current  channel  has  been  assumed  to  be  uniform  in  space.  The  physical 
p.irri  :it  sheets  haw  finite  widths  (Lc)  and  the  introduction  of  this  new  scale 
length  can  be  expected  to  make  the  results  a  function  of  (Lc/Lg).  1'or  large 
eiioih'l  Lc ,  one  can  expect  the  present  results  to  emerge.  On  the  other  hand,  for 
small  Lr,  the  local  theory  will  be  recovered.  The  precise  variation  of  \  with 
U-t  /Lg )  would  be  a:,  important  study.  We  have  carried  out  a  similar  study 
elsewhere*-  tor  the  current  driven  ion  cyclotron  mode. 

The  method  nt  solution  employed  here  was  based  on  using  a  local  expansion  of 
the  s'  function.  Direct  numerical  integration  of  the  equation  should  be  carried 
out  (the  so  called  shooting  code  method),  to  ascertain  the  accuracy  of  present 
results.  Even  within  the  analytical  method,  various  self  consistency  checks 
indicate  that  the  full  Eq .  (12)  must  be  used  in  the  strong  current  regime  (i.e. 


n 


f  :  :ui  1  1  v  ,  t  !;<  pro, i  h 1  1  i  t  y  of  t  r . :  vl1  1  1  i  r,  r  p-.Jc  s  (  i  r.  c< t  ran  t  to  ti.<- 


Iii.r::  i  ! 


c  s  t  i:d  i  t  >1  herd  should  ho  exar.i  :,ed,  takinc  Into  account  the  practical 
,'cv  i:.i  tried  f  entures  «•!'  tin  physical  dor.ai  n  ot  interest. 

Ut  had  to  lir.it  the  scope  of  this  work  due  to  limitation  c.f  resources  and  w< 
hop"  to  pursue  these  and  other  related  questions  as  and  when  further  support 
b*  cotr.i  s  acui  labl*-.  An  intcrir.  account  of  part  of  this  work  was  pr*  sentoci^'  at  the 
\u~j  Anm.ial  X.eitinp  c :  the  llasita  Physics  Li  vision  of  the  American  Physical 
S  n  c  i  e  t  v  . 
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